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Abstract — We study the two-user multiple-input single-output 
(MISO) Gaussian interference channel where the transmitters 
have perfect channel state information and employ single-stream 
beamforming. We assume that the receivers have multiuser 
decoding capability, so that when the interfering signal is strong 
enough, it can be decoded, treating the desired signal as noise, and 
subtracted from the received signal, before the desired signal is 
decoded. For this setup, the achievable rate region is a superset of 
the corresponding region for single-user decoding receivers. We 
propose efficient methods to compute the Pareto-optimal rate 
points and corresponding beamforming vector pairs, by solving 
the optimization problem of maximizing the rate of one link 
given the rate of the other link. We do so by splitting the 
original problem into four subproblems corresponding to the 
combinations of the receivers' decoding strategies - either decode 
the interference or treat it as additive noise. We utilize recently 
proposed parameterizations of the optimal beamforming vectors 
to equivalently reformulate each subproblem as a quasi-concave 
problem, which we solve very efficiently either analytically or 
via single-variable numerical optimization. We use the proposed 
methods to illustrate the effect of the strength and spatial 
correlation of the channels on the shape of the rate region. 

Index Terms — Beamforming, interference channel, interference 
cancellation, multiple-input single-output (MISO), multiuser de- 
coding, Pareto optimality, Pareto boundary, rate region 

I. Introduction 

We study a wireless system where two adjacent transmitter 
(TX) - receiver (RX) pairs, or links, operate simultaneously 
in the same frequency band and interfere with each other. 
Each TX employs tit > 1 antennas, whereas each RX 
is equipped with a single antenna. Hence, this system is 
modeled as the so-called multiple-input single-output (MISO) 
interference channel (IC) 0. We assume that the TXs have 
perfect knowledge of the local channels to both RXs and 
use scalar Gaussian codes followed by single-stream beam- 
forming. Also, we assume that the RXs have the capability 
of performing multiuser decoding (MUD) [4]. That is, when 
the interfering signal is strong enough, a RX can decode it 
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and subtract it from the received signal before decoding the 
desired signal. The decoding is done independently, since the 
RXs are located apart and there is no coordination amongst 
them. MUD capability is an important assumption because in 
principle it leads to higher achievable rates than in the single- 
user decoding (SUD) case where the RXs treat the interference 
as noise. Note that the MUD rate region is a superset of the 
SUD one, since the latter is a special case of the former. The 
achievable rate region under these assumptions is defined by 
the so-called Pareto boundary, which consists of the points 
where the rate of one link cannot increase without decreasing 
the rate of the other. 

The objective of this paper is to propose computationally 
efficient methods for finding, in a centralized way, the Pareto- 
optimal (PO) pairs of beamforming vectors which yield op- 
erating points on this Pareto boundary. These methods are 
important because they enable a fast computation of the rate 
region. Hence, they can be used to illustrate how different 
channel realizations affect its shape or in the context of large- 
scale simulation studies of interference networks. Moreover, 
they provide valuable insight on the beamforming design and 
inspiration for practical implementations. 

The capacity region of the IC is in general unknown. 
However, we know that for strong interference it is optimal to 
first cancel out the interference and then decode the desired 
signal [5 1, whereas for weak interference it is optimal to treat it 
as additive noise [6 |. Various achievable rate regions have been 
proposed and studied, e.g. in |7|-[12|. Recently, the achievable 
rate region for the MISO IC in the SUD case has been the sub- 
ject of intense studies, e.g. in fTl |-| 18 1. In 1 15], it was shown 
that single-stream beamforming is optimal for Gaussian codes. 
In ifTTI . a parameterization of the PO beamforming vectors was 
proposed based on the properties that they use all available 
transmit power and lie in the subspace spanned by the local 
channels. Alternative parameterizations were proposed in f\3\ 
and |fl2l , using the concepts of virtual signal-to-interference- 
plus-noise ratio (SINR) and gain regions, respectively. In [19|, 
the rate region for the related scenario of cooperative multicell 
precoding was characterized. 

The MISO IC with MUD-capable RXs was investigated 
in ll20l where the gain potential of MUD over SUD was 
illustrated in terms of average rate at a Nash equilibrium. 
Later, in J9], the achievable rate region for MUD-capable RXs 
was formalized and in [ 1 1 a parameterization for the PO 
beamforming vectors was proposed, extending the respective 
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result of fin . A related work appears in 1 8 1 , where a simplified 
version of the Han-Kobayashi region Q was studied and 
semidefinite relaxation was used to propose power control 
schemes and compute the corresponding rate region. 

The parameterizations fiOll . ifTTIl are useful analytical tools 
and they enable a better intuitive understanding of the prop- 
erties of the PO beamforming vectors. For SUD RXs, the PO 
beamforming vectors are obtained by trading off between the 
conflicting objectives of maximizing the desired signal power 
and minimizing the interference ifTTI . For MUD RXs, the 
trade-off is between maximization of the desired signal power 
and maximization of the interference, to enable interference 
cancellation ifTUll . Moreover, another merit of these parameter- 
izations is that they substantially decrease the dimension of the 
search space for a PO beamforming vector, from tit complex 
variables to one or two nonnegative real variables for SUD 
or MUD, respectively. However, besides the dimensionality 
reduction, these parameterizations do not directly provide a 
method for efficient computation of the Pareto boundary. The 
reason is that they only constitute necessary conditions that 
the beamforming vector of each TX has to separately fulfill, 
whereas it is pairs of beamforming vectors that yield PO 
operating points. The state-of-the-art use of the parameteriza- 
tions has been to sample the parameters, consider all possible 
combinations to generate a large number of achievable rate 
pairs, and perform a brute-force search amongst them to find 
the ones comprising the Pareto boundary, e.g., see IfTTI . 

It is desirable to devise a method which directly and 
efficiently computes PO points. Joint optimization of the 
beamforming vectors is required for this purpose. As shown 
in Ifl6l . the problem of jointly maximizing a common utility 
function and of finding PO points is NP-hard in general. 
Nevertheless, several methods have been recently proposed, 
e.g. fl4l . EIll . that apply successive convex optimization 
techniques on the vector space of the beamformers to find 
the Pareto boundary in the SUD case. The methods proposed 
in this paper achieve much higher computational efficiency by 
optimizing instead on the parameter space which characterizes 
the PO beamforming vectors of the MUD region. For SUD, 
this was attempted in 11221 . where monotonic optimization was 
utilized to find specific PO points, e.g., maximum weighted 
sum-rate points. However, the proposed method therein was 
faster than an unstructured brute-force search, but far less 
efficient than the methods that we propose herein. 

A. Contributions and Organization 

In Sec. H2 we give the system model. In Sec. [TTTJ we 
define the MUD achievable rate region and formulate the 
optimization problem that yields an arbitrary point on the 
Pareto boundary, as a maximization of the rate of one link 
given the rate achieved by the other link. In Sec. IIV44VII 
we propose very efficient methods to solve this problem, 
by combining, unifying, and improving the preliminary ap- 
proaches that appeared in our conference contributions fij, 
0. The common denominator of these methods is to exploit 
the parameterizations and equivalently recast the maximization 
problem so that it can be solved analytically or via single- 



variable optimization. The proposed method is applied on the 
following regions, whose union constitutes the MUD region: 

1) In Sec. IIV1 we propose two methods to find the re- 
gion where both RXs treat the interference as noise. In 
the first method, we equivalently reformulate the originally 
non-concave rate maximization problem as a quasi-concave 
problem of a single real variable, which we solve optimally 
with a gradient search approach. We find the Pareto boundary 
by repeating this optimization for various choices of the input 
rate. This method is novel and improves the one in ETTl in 
two ways: a) we reduce the feasible set significantly from 
the set of beamforming vectors to the parameter set defined 
in IfTTI and b) we solve a single quasi-concave optimization 
problem instead of a sequence of convex feasibility problems. 
In the second method, we use the KKT conditions to derive a 
closed-form equation which couples the parameters of the TXs 
yielding a PO pair. We find the Pareto boundary by repeatedly 
solving the cubic equation resulting from various choices for 
one of the parameters. This method was presented in fil; 
herein we give in addition a formal proof of global optimalityU 

2) In Sec. |VJ we find the two regions where one RX decodes 
the interference, before decoding the desired signal, while the 
other treats the interference as noise. We use the parame- 
terization of f£| to equivalently recast the rate maximization 
problem as a quasi-concave problem of three real variables 
and determine the solution in closed-form. This method was 
presented in 0; herein, we elaborate the proofs. 

3) In Sec. [VI] we find the region where both RXs decode 
the interference. We use the parameterization of [9| and split 
the rate maximization problem in two quasi-concave scalar 
subproblems. This method improves the corresponding one in 
in two ways: a) the number of variables is decreased from 
four real variables to a single one and b) a single instance of 
the two quasi-concave subproblems needs to be solved instead 
of a sequence of convex feasibility problems. 

In Sec. IVIII we provide illustrations of the rate regions 
for different channel properties. In Sec. IVIII1 we make some 
concluding remarks. 

The source code that implements the proposed methods 
and generates the illustrations in Sec. IVIII is available at 
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B. Notation 

Boldface lowercase letters, e.g. x, denote column vectors. 
{■} H denotes the Hermitian (complex conjugate) transpose of 
a vector. The Euclidean norm of a vector x is denoted ||a;||. 
By x ~ CJ\f(0,(j 2 ) we say that x is a zero-mean complex- 
symmetric Gaussian random variable with variance a 2 . We 
denote the orthogonal projection onto the space spanned by 
the vector x as Tl x = xx H / \\x\\ 2 . The orthogonal projection 
onto the orthogonal complement of x is JT^ = I U x , where 
I is the identity matrix. Note that for a vector y, we have 

U^y .We let f'(x) and f"(x) denote 



|y|| 2 = l|n.y|| 2 - 



'An equivalent result from a game-theoretic viewpoint was independently 
derived in 1171 and latter extended in [18 |, but without providing a rigorous 
proof of global optimality. 
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the first and second derivatives, respectively, of a function 
f(x). We define [x]% = max{x, min{x, x}}. 

II. System Model 

We assume that the transmissions consist of scalar coding 
followed by single-stream (rank-1) beamforming and that all 
propagation channels are frequency-flat. The matched-filtered 
symbol-sampled complex baseband signal received by RX; is 
then modeled as 



The achievable rate region: The rate region for the MISO 
IC with MUD capability is obtained by the union of the regions 
corresponding to all decoding scenarios, i.e., 1Z = lZ nn UlZ dn U 
7Z nd U lZ dd . Next, for each decoding scenario and a given pair 
of beamforming vectors (101,102), we define the maximum 
achievable rates QO), ifTTI . 

lZ nn - Both RXs treat the interference as noise: When 
both RXs treat the interference as noise, the maximum achiev- 
able rates of the links are [20] 



I), 



(1) 

E C" T is the 



where i,j g {1,2} and j 7^ i. In (|T), hji 
(conjugated) channel vector for the link TXj — > RXi. We 
assume that TX^ perfectly knows the direct and crosstalk 
channels, ha and hij, respectively, and that these are neither 
co-linear nor orthogonal. Also, Wi G C™ T is the beamforming 
vector employed by TX^, s,; ~ CAf(0, 1) is the transmitted 
symbol of TX^, and ~ CAf(0, of) models the receiver noise 
at RXi. The TXs have power constraints that we, without 
loss of generality, set to 1 and define the set of feasible 
beamforming vectors as W = {w € C" T | ||ti>|| 2 < 1}. The 
achievable rate for RX; depends on the received powers 



R? n (w u w 2 ) = log 2 ( 1 
R2 n ( Wl ,w 2 ) = log 2 




and 



(8) 



TZ dn - RXi decodes the interference, RX 2 treats it 
as additive noise: Since RXi decodes and subtracts the 
interference caused by TX 2 , it experiences an interference-free 
signal and the maximum achievable rate for link 1 is 



i2T(toi) = log 3 (l+Pi(«i)H 



(9) 



Pi(Wi) 



\h"w,\ 2 



and qi(wj) 



(2) 



RXi is able to decode the interference from TX 2 , considering 
its own signal as noise, if the rate of link 2 is upper bounded 
by l°g2 (l + <li( w 2)/(pi(wi) + of)) . Since RX 2 does not 
decode the interference, the rate of link 2 is also upper 
bounded by log 2 (l + p 2 (w 2 ) / (q 2 (wi) + erf)) ■ Hence the 
maximum achievable rate for link 2 is given by 



over the direct and crosstalk channel, respectively. 

In order to simplify the subsequent notation, we define the 
following channel-dependent constants. First, we define = R 2 n {w\, w 2 ) = log 2 ( 1 + min 



\hij\\ and Ki 



h'[h, 



I fay 1 1 ||fa«||)> 3 7^ *■ The latter can 
be interpreted as the cosine of the angle between the direct and 
the crosstalk channel, so when K\ = 1 or K\ = the channels 
are parallel or orthogonal, respectively. Second, using these 
constants, we define 



qi(w 2 ) 



P2{W 2 ) 



oii = ||nh y fajj|| = guKi, j ^ i, 



(3) 



n^fa« 



Pi = \\U hu hij 



9ii ~ a i= 9u V 1 - «i , 3 ^ i, (4) 
= QijKi, j 7^ h (5) 



n h«fa« 



9ii-P( = 9iidl-Ki, 3+i- (6) 



III. Achievable Rate Region of MUD Capable RXs 

In this section, we give, for completeness, the definition of 
the achievable rate region for the described scenario flOl . ifTTIl . 
We also denote the core optimization problem than we need 
to solve to find a point on the Pareto boundary. 

Each pair of beamforming vectors (wi,w 2 ) and combi- 
nation of decoding strategies (decode the interference {d) or 
treat it as noise (n)) is associated with a pair of maximum 
achievable rates. We denote R^ v (wi,w 2 ) the rate of link 
i = 1,2, in bits per channel use (bpcu), where x and y 
stand for the decoding strategies (n or d) of RXi and RX 2 , 
respectively. For each decoding scenario, i.e., combination of 
decoding strategies, we obtain a rate region by taking the union 
over all feasible beamforming vectors, i.e. 

(J {Rf{w 1 ,w 2 lRl v (w 1 ,w 2 )). (7) 



Pi{wi) + a\ ' q 2 (w{) + a\ 
(10) 

where we have used the fact that the logarithm is a 
monotonously increasing function. 

For link 2, we note that the rate is not necessarily selected to 
fully utilize the signal-to-interference-plus-noise (SINR) ratio 
at RX 2 . Actually, link 2 might hold back on its rate to enable 
RXi to decode the interfering sig nal@ 

lZ nd - RX 2 decodes the interference, RXi treats it 
as additive noise: This case is identical to 7Z dn , but with 
interchanged indices. 

K dd - Both RXs decode the interference: Both RXs 
decode the interference before decoding their desired signals. 
Since RXi decodes the interference from TX 2 , the rate of 
link 1 is upper bounded by log 2 (l + pi(wi)/of). RX 2 can 
decode the interference caused by TXi if the rate of link 1 is 
upper bounded by log 2 (l + q 2 (wi)/(p 2 (w 2 ) + erf)) . Then, 
the maximum achievable rate for link 1 is 

'P1O1) 32 (wi) 



R dd {wi,w 2 ) = log 2 I 1 + min 



°l 'p 2 (w 2 )+aj 



(11) 



By symmetry, the maximum achievable rate of link 2 is 

' p 2 {w 2 ) qi(w 2 ) 
a\ 'pi(tOi) + CT? 



R dd { Wll w 2 )=\og 2 1 



(12) 

The problem of interest is to find the so-called Pareto 
boundary of the region 7Z, which consists of PO rate pairs. 

2 This fact was not exploited in [20 Prop. 6 a)], so the description there led 
to an over-restrictive condition, hence, to a smaller achievable rate region. 
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Definition 1. A point (R*,!^) € TZ is (weakly) Pareto- 
optimal if there is no other point (Ri, R 2 ) & TZ with R\ > R\ 
and R 2 > R%. 

Graphically, the Pareto boundary is the north-east boundary 
of the region and due to Def. Q] it also includes the horizontal 
and vertical segments. In order to find the Pareto boundary 
of 71, we first find the Pareto boundaries of 7l nn , TZ dn , TZ nd , 
and TZ dd . Second, we consider as boundary of TZ the boundary 
of the union of the four 7Z xy regions. We denote by B xy the 
boundary of TZ xy and by B the boundary of TZ. 

We observe that we can find a point (R\, R 2 ) on B xy when 
the rate of one communication link, e.g. R\, is known [23 
Prop. 6.2]. The other rate, R 2 , is the maximum one that we 
can simultaneously achieve, e.g., see the dashed lines in Fig. [3] 
for TZ nn , and we find it by the following rate optimization 
problem 



maximize R 2 y (wi,w 2 ) 

(w 1 ,w 2 )e.w 2 

subject to R^ y (wi,w 2 ) = H*. 



(13) 
(14) 



Note that (Tl3T> — (fT~4l > accepts as input the coordinate R* of the 
sought PO rate pair and yields as optimal value the other 
coordinate R 2 and as optimal solution (w^w^) the enabling 
PO pair of beamforming vectors. The choice of using R* as 
input to the optimization is arbitrary. By the symmetry of 
the problem, we can choose R 2 as input and have R* as the 
optimal value. In the next sections, we derive efficient methods 
for solving (fT3l-(fT4l for all MUD-constituent regions. Since 
the logarithm is monotonic, we use the equivalent reformula- 
tion of (TT3l-(fT4l> as an SINR optimization, where the input 
parameter is 7*, i.e., the SINR (or the SNR after interference 
cancellation) required to achieve R*. 

IV. Both RXs Treat the Interference as Noise 

In this section, we compute the boundary B nn . We let 
Ri denote the maximum rate of link 1, achieved when TXi 
operates "selfishly" by using its maximum-ratio (MR) trans- 
mit beamforming vector = arg max^ e w pi(w\) = 
hxi/\\hix\\ and TX2 operates "altruistically" by using its 
zero-forcing (ZF) transmit beamforming vector w 2 F = 



P2{W 2 ) 



[20|. 



argmax W2( zw 

<?l(u>2)=0 

This combination of transmit strategies yields the PO point 

( nnn/...MR ^.,ZF\ Emn/„,,]V~~ 

\R\ V w 1 ,w 2 ),R 2 (w 1 




The rate in ( TTol l is the lowest strongly PO rate of link 
2 ll20l . Interchanging the indices in <fT3T> — (TT~6b we get the 
point (R™ 1 ,R 2 ). As illustrated by the example in Fig. [3] 
these points split B nn in three segments. The weakly 
PO horizontal (vertical) segment [(0,T% n ),(Jjg n ,R% n )] 



\[(R{ ,0),(i?i ,RT)]) is achieved when TXi (TX 2 ) uses 
the MR beamforming vector and TX2 (TXi) uses the ZF 
beamforming vector, adapting the transmit power in [0, 1]. 

The remainder of this section focuses on the strongly 
PO segment between (i?™ i? 2 ) and (#1 iMa)- Inserting 
© into (fT3l-(fl4l and equivalently reformulating the rate 
maximization to SINR maximization, we obtain 



maximize 
subject to 



P2{w 2 ) 
q 2 (w 1 ) + o\ 



It- 



(17) 



(18) 



qi(w 2 ) + of 

From fl5i , we see that constraint ( fl8l ). hence the optimization, 
is feasible when 7^ < g\ l ja\. The formulation (fT7b-(fT8b is 
non-convex since the objective function ( fTTb and the equality 
constraint ( fT8l consist of fractions of quadratics. 

In Sees. IIV-AI and IIV-B1 we propose two methods that 
find very efficiently the global optimal solution. The main 
difference between these methods is the input required to yield 
the entire boundary; in the first method it is different choices 
for one of the PO SINR values, whereas in the second method 
it is different choices for one of the PO beamforming vectors. 
Both methods have computational complexity that is constant 
in the number of transmit antennas. The method in Sec. IIV-AI 
can be interpreted as solving an underlay cognitive radio 
problem where the secondary user, here link 2, maximizes 
its rate under various quality-of-service constraints for the 
primary user, here link 1. The interpretation of the method 
in IIV-BI is that TXi fixes its beamforming strategy and TX2 
seeks the best-response strategy to end up at a PO point. 

A. Numerical Method 

The numerical method proposed in this section is a two- 
fold improvement of the one we presented in ETI . First, 
we exploit the parameterization [11] of the PO beamforming 
vectors to reduce, without loss of optimality, the feasible set 
from W 2 , i.e. a bounded convex set in M 4llT , to the bounded 
positive quadrant. Second, we equivalently reformulate the 
optimization problem to further reduce the feasible set to a 
line segment in R. The resulting problem is a scalar quasi- 
concave problem that needs to be solved once to yield a PO 
point, whereas the parameterized convex formulation in 12D 
required several iterations of the bisection method. 

From 1 11 Corollary 1], we know that the PO beamforming 
vectors of the 7Z nn region can be parameterized as 



w l (x i ) 



l-x 



(19) 



where < Xi < 1 for i,j — 1,2 and j 7^ i. Note 
that ( [19] l consists of a single nonnegative real parameter per 
beamforming vector. Inserting ( fr9l into (O, we get 



Pi{Wi) 



I ha j 



= CliXi +Cti\/1 



(20) 
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|n hi;j hu\ 



2 2 
9ij x i ■ 



(21) 



From ifTTl . we know that, for PO points, Xi < cui/gu = /tj. 
This value maximizes ( 120b and corresponds to the MR transmit 
beamforming vector. Further increase of Xi will just increase 
the interference and decrease the desired signal power. 

Inserting ( f20b and (fJTJ into (TTTb — (TT~8T> and performing 
straightforward algebraic manipulations, including taking the 
square root, we get the equivalent reformulation 



maximize 

Q<X\ ,X2 < 1 



subject to 



Ot2X2 + OLl\J\ — X 



(22) 



(23) 



where the constants a%, 67%, a.2, cb) <?12, <?2i are all positive, as 
defined in Sec. [TT] We further simplify the notation defining 



Ui(Xi) = OL i X l + OliJl 



h3 



i = 1,2 and 
1,2, j^i. 



(24) 
(25) 



It is straightforward to verify that Ui(xi) is concave and non- 
decreasing for X{ < Ki- Moreover, since Vi(xj) is a norm, it is 
a convex and non-decreasing function in X{. Then, we make 
the following observation: 



Lemma 1. The objective function (I22l i and the left-hand- side 
(LHSj of the constraint (I23l l are quasi-concave functions. 

Proof: Note that Ui(xi)/vi(xj) > c is equivalent to 
cVi(xj)—Ui(xi) < since Vi(xi) > 0. Since cvi(xj) is convex 
for all c > and Ui{xi) is concave, cvi(xj) — Ui(xi) < 
defines a convex set. Hence, we can conclude that the objec- 
tive d22l and the LHS of constraint ( 1231 ) are quasi-concave 
functions [24, Ch. 3]. ■ 



Due to the equality in (1231 . the problem (I22l-(l23l is not 
quasi-concave as it stands [25 1, but in the following we 
equivalently reformulate it into a quasi-concave problem in 
one scalar variable. We solve equation ( f23l for x%, keeping 
the positive root, as 



X2 



fl2l7l 



(26) 



Since a function of the form y/t 2 — a is concave and non- 
decreasing for t > yfa, a > and ui(xi) is concave and non- 
decreasing for xx < Ki, we conclude that w(x%) is a concave 
and non-decreasing function of x% < k,\ ll24l Ch. 3.2]. 

The fact that < X2 < ^2 introduces lower and upper 
bounds on x\. Since X2 > 0, it follows from (f26b that 
Mi(^i) > \fli°~i- Then, from (l24l and X\ > 0, it follows 
that we must have xi > x lt where 



S.i =max^0, -j-J^l- -r\/9u- lief 
9u 9u 



(27) 



Note that x_ x is real for 7* < g^/0 2 . Furthermore, for the 
upper limits we must have X2 — w(xi) < K2 and x\ < m, 



which imply x% < X\, where 



Oil 

9li 



V^Tl(«2ff21 +°"l) 



X\ = < 



Oil 



9 2 u 



(28) 



where 7 f R = g\ x l{glx^\ + °\) is the SINR of link 1 when 
both TXs use the MR beamforming vector. 

Inserting d26l i in (l22l . along with the lower and upper 
bounds (|27T > and (|28T >. respectively, gives us the unconstrained 
scalar optimization problem 



maximize 



ita(tw(xi)) A 



s(ii). 



(29) 



U2OB1) 

Note that the objective function corresponds to the square root 
of the SINR of link 2, i.e, s(xi) = y/yz. Next, we study its 
properties and prove that it is quasi-concave. 

Lemma 2. A stationary point of s{x\), if it exists, must be a 
maximum. Hence, the function must be quasi-concave function. 

Proof: We know that w{x\) is concave and non- 
decreasing for x\ < Ki, V2(xi) is convex and non-decreasing 
and 1*2(^2) is concave0 At a stationary point x\, the first 
derivative of s(x\) satisfies s'(x*) = 0, where 

w' {x\)u' 2 {w{x\))v 2 {x\) - U2{w{x\))v' 2 (x\) 



s'(x\) 



v 2 {x\) 



(30) 



(31) 



We observe that for a stationary point, we have 

u , {w{<)) = > 

W'{XI)V2(X*) 

since u 2 (w(x*)), v' 2 (x*), w'(x*), and V2 (x*) all are non- 
negative for x* < K\. We have equality in (fJTJ) only if x* = 0, 
but this case is treated separately. Now, it is straightforward 
to see from the curvatures of Ua^a), ^2(^1), and w(xi) and 
( l30l l that a stationary point satisfies s"(x*) < 0, where 

= -2^-Aw"{x\)u' 2 {w{x\))v 2 {x\) + 
v 2 \ X l) 

+ ( w ' ( x i)) 2u 2( w ( x i)) v 2(xl) - u 2 (w(x*))v 2 (x*)), (32) 

since all three terms are non-positive. The term 
U2(w(x*))v' 2 (x*) = only if x\ — 0. Hence, if there 
exists a stationary point, it must be a local maximum. Since 
s(x\) is continuous differentiable for x\ > x x , that local 
maximum must be the global maximum. Hence, we have 
shown that s{x\) is quasi-concave and it has no saddle points. 

■ 

Since problem d29l has a single real variable and the 
objective function is quasi-concave, the optimum solution can 
be found very efficiently. Since there are no saddle points, a 
steepest gradient method can be used. In Tab. [I] we propose 
an algorithm that computes strongly PO points of TZ nn . As 
input, the algorithm requires the channel vectors and noise 

3 We could have used the fact the U2(w(xi)) is also non-decreasing for 
xi < xi, to obtain a simpler proof. However, in Sec IVII we need this more 
general case. 
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variances, and the number M of requested boundary points. 
The output is stored in the vectors ri,r 2 <E M. M . The rates 
in r\ are obtained by uniform sampling over the interval 
In line 4, we compute the end point #2™)- 
For each boundary point, we compute the lower and upper 
bounds x_i and x\, respectively, and ensure that the solution 
lies in the interval [ ]. In line 9, the x* corresponding to 
the previously computed point is used as initial value for the 
next point on the boundary. The reason is that we expect that 
the solution will not change significantly for two nearby points. 
In lines 1 1— 16, we find the optimal solution x\ by a gradient 
ascend method. In each repetition, we compute the derivative 
s'(x\) and then find a step length t, by backtracking line search 
ll24l Ch. 9.2]. This is repeated until the improvement from the 
previous iteration is smaller than some predefined tolerance 
e. Since s(xi) is quasi-concave, its derivative can be small 
without being close to the optimum. Hence, e has to be chosen 
very small. In line 23, the end point (R 1 ,R^ n ) is computed. 



Input: hij, af, i,j = 1,2, M, and e 
Output: B nn given by vectors r±,r2 S K M 
ri = \R1 n : (R" n - R™)/(M - 1) : 7?"" 

r-a(l) = RT< x *i=0 
for k = 2 : M - 1 

7 * = 2* , i('=) - 1 

Compute x_i using J27t 

Compute x\ using j28t 



, =[x[ l> +ts'(x^)]%i 



x\ = \x 
I = 
repeat 

Compute s'(x^ ) 

Determine step size t 

,.('+!) _ r_W 
x 1 — 

li- 1+ 1 
until \s{x[ l) ) - < e 

x* = 

Compute xi = w(x^) using <26t 
Compute id* = Wi(x*) using H9\ 
Compute i?2 = K^*^ 10 *' w 2) usine {8j 
r 2 (fc) = RJ 



end 

r 2 (M) = R% n , x* 



TABLE I 

Numerical method to compute B nn 



B. Closed-Form Parameterization 

In this section, we use the Karush-Kuhn-Tucker (KKT) 
conditions of the optimization problem (f22l>— (l23l> in order to 
derive a closed-form relation between the parameters of the 
beamforming vectors that jointly yield a PO rate point. A pre- 
liminary version of this method was presented in (TJ; herein, 
we elaborate the derivations and provide a proof of global 
optimality. The latter is achieved using the parameterization 
(fl9l l. whereas a different parameterization was used in |Q]]. 

In general, the KKT conditions only provide necessary 
conditions for global optimality. However, we show that for 
this specific problem, the KKT conditions are also sufficient. 
Towards this direction, we relax the equality constraint (1231 to 
a lower-bound inequality. As evidenced by the dashed lines in 
Fig. [3] this relaxation is tight at optimum; hence the relaxed 



problem is equivalent to the original one. Then, due to Lem. Q] 
the relaxed optimization problem (t22b — (f23l> falls into the class 
of quasi-concave problems 11251 . Th. 1 in (25 1 gives a number 
of sufficient conditions for global optimality of the solution to 
the KKT conditions of a constrained quasi-concave program. 
It is enough that one of these conditions is satisfied. Condition 
a) is that the gradient of the objective function should have 
at least one negative component for a solution that satisfies 
the KKT conditions. By simple inspection of the objective 
function (|22"||. it follows that: 



Lemma 3. The objective function ( 1221 i is decreasing with X\ > 
0, for fixed X2- 

Hence, due to Lem. [3] the relaxed version of the problem 
(F22l -(l23l satisfies condition a) of Th. 1 in ll25l . Then, from 
Lem. Q] Lem. |3] and 11251 Theorem 1], we have the following 
result: 

Proposition 1. The KKT conditions of the relaxed problem 
(122l l- (1231 > are sufficient conditions for global optimality. 

For notational convenience, we make the bounding con- 
straints on xi implicit, i.e., we declare a solution of the KKT 
conditions feasible only if it adheres to the bounding con- 
straints. The case of = corresponds to ZF beamforming, 
which was treated in ( TTol ). The Lagrange function of the 
relaxed (f22b — (f23T> is 



C(xt,X2,p) 



U2(X 2 ) 



/' 



ui(ari) 



'7i 



(33) 



v 2 {xi) \vi(x 2 

where the Lagrange multiplier p is non-negative. Hence, the 
KKT conditions become [24| 



dC 

dx\ 
dC 
dx 2 



V Vi (X2 ) 
V' 2 {XI)U2{X2) 

vl{xi) 
u' 2 (x 2 ) 



'7i 
u[(xi 



v[{x 2 )ui(xi) 



= 0, 
= 0, 
= 0. 



(34) 
(35) 
(36) 



V2{xi) ^ Vf(x 2 ) 

In (T34ll-(l36ll. we avoided explicitly including the primal fea- 
sibility constraints, since (t22b — (f23l> is always feasible if we 
choose 7^ < g\ijo\. We use the KKT conditions to find a 
relation between the parameters x% and X2 of the beamforming 
vectors that jointly yield PO points. Since we are not looking 
for a specific PO point, we can disregard condition d34l l. As 
long as the other conditions are satisfied, we can find the 7* 
that satisfies 041 . 

When x\ > 0, we can verify from d35b that p > 0. Then, 
we use (1331 and (f36b to solve for /1. Equating the solutions, 
we get the relation 

U 2 {X2)vl(X2) V , 2 (Xi)v 1 (X2)u2(X2) 



V , 1 {X2)v 2 (Xi)u 1 (x 1 ) 



U[(X!)V% (Xi) 



(37) 



By collecting all functions of x\ and X2 in the LHS and RHS, 
respectively, we equivalently rewrite (T3~7l as 

,/ \ A U[(xi)v%(xi) _ V' 1 (X2)V1(X2)U2(X 2 ) a / s 

J[ x l) — I, \ 1 \ I \ — 1, \ 21 \ ~ 9( x 2) 



v' 2 (xi)v 2 (xi)ui(xi 



u' 2 (x 2 )vl (x 2 ) 



(38) 
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The LHS and RHS of (I38t are functions of only x\ and x 2 , re- 
spectively, which we denote as f(x\) and g{x 2 ), respectively. 
In order to find a PO point, we fix ii at a specific value x * 
and then solve g(x 2 ) — f{x\) to get xj. 

Due to the square roots, it is complicated to solve for x 2 
as it stands. Instead, we use the alternative parameterization, 
that the PO beamforming vectors are linear combinations of 
the MR and ZF beamforming vectors ifTTl Corollary 2], i.e., 

X iW f R + (1 - \i)wf 



PO 



(A, 



X t wf R + {1-Xi)w 



ZFI 



(39) 



where Xi £ [0, 1]. To go from the parameterization in (fT9l l to 



that in 



we use the mapping 

OLi 



Xi = (f>i(Xi) = 



A, 



9» \\X t wf R + (1 ~ Xi)w 
Xi 



ZFI 



(40) 



where pi = 1 — yl — k 
Pi < 1 and 

d4>i _ on 
dX 



9u v / 2p i X'f - 2 Pl Xi + 1 

2 Since < Ks < 1, we have < 



1 - PiXi 



> 0. 



(41) 



<to (2pA 2 - 2 Pi Xi + l) 3 / 2 

Hence, it follows that (l40l > is a one-to-one mapping between 
Xj and Xi. Therefore, it follows that the problem of solving 
5(^2) = f(xi) with respect to x 2 is equivalent to that of solv- 
ing <?(</> 2 (A 2 )) — /(0i(A|)) with respect to A2. By inserting 
(l40l into OSt . we equivalently write g((f> 2 (X 2 )) = f((f>i(XD) 
as 

A 2 (l - ,o 2 A 2 )(p 2 A 2 + (1 - p 2 )) 



(1-A 2 )(p 2 (2- P2 + 2C 2 )A 2 



2P2Q2X2 + C2) 

(42) 

i,j = l,2, j 7^ i. In (l42l . we see that 
<?(</>2(A 2 )) is a fraction of cubic polynomials. Since /(<t5i(A{)) 
is a constant, we write d42l as the cubic equation 



where = a 2 /g 2 . 



Pi 



(43) 



(44) 



C3A2 + c 2 A2 + ciA 2 + c = 0. 
The coefficients of the cubic equation d43l are 

co = -C 2 /(Ai), 

c 1 ±(l + 2p 2 )( 2 f(X* 1 ) + (l-p 2 ), 

C2 = ~P2{2 - P 2+ ^P2C2) f(MK)) + 

c 3 ^p 2 (2-p2 + 2C2)f(MK))~Pl 
Cubic equations can be solved in closed form [26]. The roots 
of ( |43T > are three candidates for A£. Since A^ £ [0, 1], we have 
the following three cases. 

. A{ = 0: From dT6b and d39b we know that Ag = 1. 

• < A{ < 1: We find the roots of ( l43l and keep the roots 
that satisfy the feasibility constraint, that < A 2 < 1. 
We can potentially have more than one feasible root, but 
due to PropJTJ we know that all feasible roots yield a 
PO solution^ 

• A{ = 1: Again, from (O and ([39]) we see that A£ = 0. 



4 In [171 and 1 18], it was not made clear whether all feasible solutions to 
the corresponding cubic equation are optimal or not. Especially, equation (25) 
in 1171 provides only necessary conditions for Pareto optimality 1231 Ch. 6]. 



The overall method to compute the entire B nn is summa- 
rized in Tab. [TT] Once we have found the coefficients in Sec. 
Ull the complexity is constant in the number of antennas. 



Input and output: same as in Tab. |I] 



for A 



[0 : 1/(M - 1) : 1] 



Compute /(</>i(AJ)) according to )38t and 1401 

= roots of cubic equation <43t that are in [0, 1] 
Compute the rate point(s) using {8), U9t . and )40t 



end 



TABLE II 

Closed-form method to compute B nn 



V. Only one RX Decodes the Interference 

In this section, we propose a closed-form method to com- 
pute B dn . A condensed description of this method was given 
in 0. 

Inserting the rate expressions (O and ( fTOb in the optimiza- 
tion problem (fT~3b — (TT4l > and equivalently reformulating the rate 
maximization to SINR maximization, we obtain 



maximize 

(iui ,w 2 )ew 2 



qi(w 2 ) 



P2{W2) 



Px{wx) + a\ ' q2{w x ) + a\ 



subject to pi{wi)/a\ = 7*. 



(45) 
(46) 



As with (ITTb— (fl~8T>, we see that d45l>— (f46l> is feasible when 7* < 
g\ijo\. The formulation (I45ll-(l4"6ll is nonconvex, because the 
objective ( f45b is the minimum of two fractions of quadratic 
functions and d46b is a quadratic equality constraint. 

In iflOl . it was shown that the PO beamforming vectors of 
the lZ dn region can be parameterized as 



wx(xx,yx) = X\ 

W 2 (x 2 ) 



X2-, 



rih 12 /iii 

|n hl2 h,n| 

\N-h 22 h2i\ 



2/1- 



1-x 2 



|n^>2i| 



(47) 
(48) 



where [x\, y\) £ Q = {(x, y)\x, y > 0, x 2 +y 2 < 1} and x 2 £ 
[0, 1]. Note that this parameterization is different from ( |T9l 
and d39l of lZ nn . An interpretation stemming from d47li-(l48Tl 
is that on the Pareto boundary TX 2 uses full power, whereas 
TXi may not. Inserting (l47b-d48T) into ©, we get 

2 



Pi(wx) 
92 Oi) 

P2{W 2 ) 

qx (w 2 ) 



(xi ||n hl2 ft,n|| 

{aixi + aij/1) 2 

2 |^12^1l| 2 



•»i n h hu 



ll n h. 12 ^iill 

2 \h2 2 h 2 i\ 2 



l|nh, 22 /i 2 i M 
x 2 ||n h22 ft 2 i 



_ 2 2 

2 ~ 9x2 x \i 



— 22 
2 ~ 922 X 2i 



(49) 
(50) 

(51) 



$2X2 + /W1 



1 -x 

2 



(52) 



where the parameters olx, otx, (3x, fix, gi2, g 2 2 are positive, as 
defined in Sec. [TTJ When x\ increases, both the power of the 
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desired signal d49b and the interference (T50b increase, whereas 
yi only increases the power of the desired signal. When x 2 
increases, the desired signal power ( BIT ) increases. 

Inserting d49b— (f52b in (|43T> — (l46b. replacing the constraint 
(l46l > in the denominator of the first fraction of (l45l ). and taking 
the square root, we equivalently obtain 



/3 2 x 2 + fe-y/l ~ x 



maximize 

(«,si)eS 



subject to aiXi + = \/J*cTx- 



922X2 



(53) 



By Lem. Q~| we see that the two fractions in the objective 
function (l53l are quasi-concave. Since the minimum of two 
quasi-concave functions is quasi-concave and (l54l ) is linear, it 
follows that (I53l-(l54b is a quasi-concave problem. 

We solve (l53l-(l54b in two steps. First we solve for (x±,y%) 
and then for x 2 . We note that x% appears only in the second 
fraction of ( f53l and in (l54l l. whereas y% appears only in ( 1541 ). 
The second fraction of d53l is monotonously decreasing with 
X\, for fixed x 2 , so we maximize it by minimizing x\, subject 
to the constraint ( 1541 . i.e., 



minimize 

Oi,yi)eg 

subject to 



yi = -—xi 



(55) 



(56) 



The solution of this problem can be found by inspection, 
noting that the feasible set is the segment of the line d56l > 
in Q. When \/ 7* of / &i < 1, this line segment does not 
intersect the unit-radius circle; hence, the optimum value is 
x\ = 0. The interpretation of this solution is that TXi uses 
the ZF beamforming vector, but contrary to the lZ nn case, 
it may not use full power in order to make interference 
cancellation possible. When ^/y\a\ja\ > 1, the line segment 
intersects the unit circle in two points; hence, the leftmost 
is the optimum. Inserting d56l l into the quadratic equation 
of the unit circle, it is straightforward to determine that 

x\ = faiVrfof - «lV9u l* a i) /Sii> where we have 
used the fact that g 2 x — a\ + a\. The interpretation of this 
solution is that TXi uses full power in this case. 

Given the optimal x*, the optimal y* is determined by 
and the problem (T53l>— (f54b only depends on x 2 , i.e., 



maximize mm 

0<x 2 <l 



P2X2 + (3 2 y/l - x 



522^2 



(57) 

In order to simplify the notation, we define the constants 

a = 922/\f{x\) 2 gl 2 + of, b = /WW (7* + 1), and c - 
$2/ y/vidi + !)■ Using these constants, we write ( 1571 ) as 



maximize 

0<rc 2 <l 



bx 2 + cy 1 — x 2 , ax 2 



(58) 



Depending on the values of a, b, and c, we get three different 
cases for the objective functions in d58l ). as depicted in Fig. Q] 
In Fig. Q] (a), we have a < b and it is clear that the optimum 
is at x 2 = 1. The interpretation is that TX2 uses the MR 



beamforming vector. The difference between the cases in Fig. 
Q] (b) and Fig. Q] (c) is whether the intersection of the curve 
with the straight line is to the left or right of the maximum 

zy/l — x\ is maximized for 



The intersection of bx 2 + cyl — x\ with 



of the curve. The curve bx 2 

x 2 = b/Vb^+ 

ax 2 is happening for x 2 = c/y^c 2 + (a — b) 2 . In Fig. [U (b), 
the intersection is to the right of curve's maximum; hence, 

b 



> 



y/c 2 + (a - b) 2 ~ Vb^T. 



^ ab <b l + c z 



(59) 



From Fig.Q](b), we see that the optimum is at the intersection. 



(54) Hence, we have x 2 = cj y/ c 2 + (a — b) 2 . For the case in Fig. 



Q] (c), we see that the optimum is at x 2 = b/y/b 2 
fi/g 2 \. For the cases depicted in Figs. Q](a) and (b), we have 
ab < b 2 + c 2 and the solution lies on the line ax 2 . Therefore, 
we have j 2 = (ax 2 ) 2 . In Fig.[T](c), we have ab > b 2 + c 2 and 
the optimum lies on the curve bx 2 + cyl — x\. Therefore, 
we have = b 2 + c 2 . The interpretation is that TX2 uses a 
beamforming vector in the direction of the crosstalk channel. 
Note that in Fig. Q] we depicted the scenario of b > c. The 
above analysis does not change if b < c. We summarize the 
solutions of (T55T>— <f56b and ( 1571 ) in the following proposition. 



(a) 



(b) 



(c) 



Pi 




•t'2 



X 2 



x 2 



Fig. 1 . Illustration of the three different cases in the proof of Prop. f3] The 
optimal solution is marked with a star, (a): a = 1, 6 = 1.2, c = 0.9, (b): 
a = 1.5, b = 1.2, c = 0.9, (c): a = 2.5, b = 1.2, c = 0.9. 



Proposition 2. The optimal solution of (I55ll-(l56ll is 

(60) 



x\ = max \ 0, -J- -J yfal - -±- <J g\ t - y\a\ 



9u 



Vi = (\J 7*0-1 - axxlj lax- 
Then, the optimal value of (145b is given as 

922( x 2) 2 

I (a 

72 



{x\9i2) + <r. 

921 



3 , a<b + c 2 /b, 



(61) 



(62) 



^(7!* + l)^ 



a > b + c 2 /b 



for 



Xo = < 



1. 



a < 6, 



(63) 



c/Vc 2 + (a - 6) 2 , b<a<b + c 2 /b, 
4 /3 2 /.92i, a >b + c 2 /b, 

The optimal (w *, w 2 ) is obtained by inserting ( 1601 ) — ( f6TT > and 
463) into 

Prop. [2] provides a scheme for finding in closed-form a 
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point on the Pare to boundary B dn , by providing 73 as an 
explicit function of 7*. In Tab. Hill we summarize the proposed 
algorithm for computing the entire B dn . 



T: Input and output: same as in Tab. |I] 

2: r x = [0 : HT/{M - 1) : i?f ] 

3: for k = 1 : M 

4: 7* = 2 7- i( fc) - 1 

5: Compute and y * us ing )60t and j6U 

6: Compute aiS using 163) 

7: Compute lUj and wt, u sing j47t and (48) 

8: Compute 7J using j62t 

9: r 2 (fc) = log 2 (l+ 7 *) 

10: end 



TABLE III 

Closed-form method to compute B dn 



Contrary to the 7Z nn case, we can obtain the weak, i.e., 
vertical and horizontal, parts of B dn by using the algorithm in 
Tab. |nl] The reason is that by using the parameterization (|47|>- 
d48l i. we can set Pi(iVi) = 0, which is not possible for lZ nn . 
Moreover, it is of interest to analyze the largest value, 73 , 
that 7J can assume, i.e., when 7J = 0, since this brings some 
insight to lZ dn . At this point we have a = g 22 /cr 2 , b = (3 2 /(Ti, 
and c = $2/^1- Therefore, for 7J = we have, 



if 



522 



= < 



M9I2 



921 



+ (522CT1 - f3 2 a 2 Y 



.922 
02 

Ih 

.922 
02 



< 



< 



§1 
.922 

c 2 



< 



.921 



> 



921 
K 2 (7l ■ 



(64) 

Since /3 2 = K 2 g 2 i, the first case of ( f64-b corresponds to the 
scenario where the crosstalk channel ft.21 is strong compared 
to the direct channel h 2 2, and the spatial correlation between 
h 2 \ and h,22 is large. For this scenario, it is optimal for TX2 
to use the MR beamforming vector. The third case of (l64l 
corresponds to the scenario where h,2i is weak compared to 
h,22, and the spatial correlation between ft.21 and ^22 is large. 
For this case, TX2 has to prioritize the SINR at RXi and it 
uses a beamforming vector in the direction of h,2i- The second 
case of ( f64-b is somewhere in between the previous extreme 
cases. For this case, TX2 chooses its beamforming vector such 
that both the receivers get the same SINR. To conclude, the 
highest rate link 2 can achieve in 7Z dn is R 2 = log 2 (l+72 )■ 



VI. Both RXs Decode the Interference 

In this section, we propose a computationally efficient 
numerical method to compute B dd which is similar in logic 
to the one given in Sec. IIV-AI for B nn and also utilizes 
intermediate results from the method in Sec. [V]for B dn . The 
method proposed herein improves the corresponding one of 
J2] in two ways: a) the number of variables is decreased from 
four real variables to a single one and b) a single instance of 
two quasi-concave subproblems needs to be solved instead of 
a sequence of convex feasibility problems. 

Inserting the rate expressions (fTTT i and (fT2l in the optimiza- 
tion problem (TT~3T> — (TT~4T> and equivalently reformulating the rate 



maximization to SINR maximization problem, we obtain 

' p 2 (w 2 ) qx(w 2 

maximize ' 
subject to 



cr 2 

Pi(t»i)/a? > 7i, 
92(^1) 



> 



7l: 



(65) 
(66) 
(67) 



P2O2) + <jI 

where d66i>-d67b follow from the so-called epigraph formu- 
lation of CD}, see [24, Ch. 3]. The formulation <fg51)-<f571) is 
nonconvex, since it has constraints that consists of fractions 
of quadratic functions. 

In 1 9], it was shown that the PO beamforming vectors of 
the lZ dd region can be parameterized as 



w i (x i ,y l ) = Xi 



(68) 



for i,j = 1,2 and j ^ i, where (xi,yi) £ Q. Note that this 
parameterization is different from ( fT9l and d39b of 7Z nn and 
d47l-d48l of K dn . Inserting d68]) into ©, we get 

\h H h 12 

p l (w i ) = x t - 



! J 



l n ^ii^j II 
/' ll n ^^ll 



2 2 
9ii x i 1 



(69) 



(70) 



where the parameters gu,l3i,j3j are defined in Sec. [TT] From 
d68l l. we see that the PO beamforming vectors of both TXs 
do not necessarily use all available power. However, without 
loss of optimality, we can assume that full power is used 
at optimum. This is so because increasing iji increases the 
interference d70l ) but does not affect the desired signal power 
d69l . The effect of increasing yi, beyond the optimal solution 
y*, is to only make the constraint ( f67b looser at optimum, 
without decreasing the objective value (f65t . The interpretation 
is that we can increase the interference arbitrarily, since it will 
be cancelled by the RXs. Hence, we can increase it until the 
power constraint is met with equality, i.e. set y. L = y/l — x\. 

Inserting (l69l-(l70li in (I65l)-(l6"7l). taking the square root 
of objective and constraints, and introducing the nonnegative 
auxiliary variable z, we equivalently obtain 



maximize 

0<a;i,a;2<l, z>0 

subject to 



922X2I02 > z, 



P2X2 + /?2\/l - x\ 

Vgli^l - 

gnxi/ai > 
PiXx+Px^Jl-xj 

\/9 l 22 X 2 + °2 



u l 
^1; 



> 



7i- 



In order to further simplify the notation, we define 



Ui{xi) = fiiXi 
Vi{xi) = yfg. 



1 



(71) 
(72) 

(73) 

(74) 

(75) 

(76) 
(77) 
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Problem d7B-CZ3 is feasible for 7? <E [0,jf d ]. The rate 
R dd = l g 2 (l + jf) is the highest rate of link TXi ^RXi 
that can be decoded by both receivers, achieved when TX2 
does not transmit. We determine jf d by solving 



maximize mm 

0<a;i<l 



511X1 P1X1 + (3iy / T^xJ 



0"! 



(72 



(78) 



Problem (F78b is similar to d58l ). so we can solve it using d64b . 
By interchanging the indices in (F78T >. we can find 73 d . 

Due to the variable z, ( l73l l does not define a convex set. 
Hence, (l7lTi-(l7"5Tl is not a concave nor a quasi-concave problem 
as it stands. But, by using the epigraph formulation, (f7TTi- 
(Tt3T > can be equivalently reformulated into a quasi-concave 
problem. By studying the KKT conditions of (l7lTi-(l7"5ll. we 
identify two cases and apply to each of them techniques 
introduced in Sees. IIV-AI and [V] respectively. The fact that 
the gradient of the Lagrange function of (|7T1i-(|75Ti vanishes at 
the optimum, gives the KKT conditions 

= 1- Ml - M2 = 0, (79) 



dz 

dC v' 1 (xi)u 2 (x2) . gn . u[(xi) 

a — - — \ 1- ^3 h Hi—, — r 

0x1 v((xi) ox vi\xi) 

dC g 2 2 . u' 2 (x 2 ) v' 2 (x2)u2(x 2 ) 

a — = A*i 1" M2-^t — r - rjT — \ - 

ox 2 <y 2 v\{x\) viy x 2) 



= 0, (80) 



0, (81) 



where Hi > 0, i € {1, 2, 3, 4} are the Lagrange multipliers of 
constraints (F72b — (F73T>. respectively. First, we observe that we 
can have /13 = L14 = only when x\ = 0. This is the case 
only when 7* = 0. Hence, for every other point we have either 
/13 > or LI4 > 0, corresponding to equality in (1741 and d75l l. 
respectively. Next, for each case, we change the corresponding 
inequality in (I7l1i-(l7"5ll to equality. We solve the two programs 
separately and compare the solutions. The solution with the 
highest optimal value will yield the optimum of (l7TTi-(l75ll. 

For the case of equality in (l74l . it immediately follows that 
x* = g\\- From (|75| ). we see that the problem is 

feasible only if /3ix\ + fixy/l - {x\) 2 > \J~f\o\. Given that 
d7Tt-d73t is feasible with equality in ( 1741 . we have the problem 



maximize 

0<a:2<l 



subject to 



g 2 2X 2 u 2 {x 2 ) 



a 2 



X 2 < 1 - T =Ju\{x\) --i\u\. 

922 V7t v 



(82) 
(83) 



The structure of d82l-d83l is similar to d57l . The only dif- 
ference is the extra constraint ( l83l which yields a tighter 
upper bound for x 2 . Hence, we can use Prop. [2] to find x\ 
by using coefficients a = 322 /02, b = ft/vT? + ^1 > an( ^ 
c = + a \ ' m pl ace °f a ,b, and c ' respectively. 

For the case of equality in ( 1751 ), we get 

x 2 = ^ 7 =J^(xt)- 7l V 2 2 ±w( Xl ). (84) 

#22^7^ v 

Inserting ( l84l in (TTTb— (T75l) yields the problem 



maximize 

0<a;i<l 

subject to 



min{si(xi),s 2 (x 1 )} (85) 
guxi/ai > y/tf, (86) 



where 



si(a;i) = g22w(xi)/a 2 , 
s 2 {xi) = u 2 (w(xi))vi{xi), 



(87) 

(88) 
(89) 



and the constraint dSTl corresponds to a; 2 < 1. Next, we 
determine the properties of (l85l-d87l). By revisiting the KKT 
conditions (T79b— (fSTb. we see that if ^4 > and x\ > K\, 
then we must have /13 > as well. This follows since 
Ui{x\) < for x\ > k±. But the case of H3 > is 
already covered by d82l-d83l. Hence, if 7* > g\ l K\/a\, 
it is enough to solve d82b— (f83l). For this case, we set the 
optimal value of (I85l-(l87l to zero. Since u\(x\) is non- 
decreasing for xi < m, the constraint d87l becomes convex. 
It is feasible when ui(0) = Pi < y/jtig^ + ^D- Assuming 
it is feasible, d87b yields an upper bound on the optimal x%; 
namely x\ <X\, where 



x\ 



9I2 



Ki, 



(922 



^r\/g 2 i2-(9h + CT l)7i*, 

#12 



7i*< 
7* > 



9 12 

2 
3l2 

922+"2 ' 



(90) 

For the lower bound on x%, we have constraint ( l86l and the 
fact that ui(xi) > 7* cr^, which yield x\ > x lt where 



a m Pi 

9h^~9h 
0, 



5?2 



7* > 



5l2 

(91) 



Next, we show that the objective function ( l85l l is quasi- 
concave for xi g [x.i,xi]. First, we note that the minimum 
of two quasi-concave functions is quasi-concave l24l Ch. 3.4]. 
The function d88l is on the same form as d26l i and hence, it 
is concave. So, if (T89T ) is quasi-concave, then d85l ) is quasi- 
concave. We note that (f89b has the same structure as the 
objective function of d29i >. Hence it follows from Lem. [2] that 
d89l is strictly quasi-concave. 

Since d85b is strictly quasi-concave, we can use a gradient 
method similar to the respective one for B nn , presented in 
Tab. U The proposed method to compute B dd is sketched in 
Tab. IIVI and differs in the following points to the one for B nn . 
First, for B dd we have to solve two optimization problems, 
which we do separately. We denote x* 1} x 2 % and 7^ the 
optimal solution and value, respectively, of d82l-d83l. and 
x* 2 , X22 and 7I2 the optimal solution and value, respectively, 
of (f85l>— (f87l>- The solution to the subproblem that yields the 
highest optimal value is declared the solution to (T7Tb— (T75b- If 
both subproblems are infeasible, we set the optimal value to 
zero. Second, we maximize the minimum of two functions. 
Therefore, in lines 14-18, we let the gradient of the objective 
function d85l ). denoted by A, at a point x\, take as value the 
minimum of the derivatives of functions si(x\) and s 2 (^2)- 
In line 22, we use s{x\) = min{si(xi), 52(^1)} • Except for 
these points, the method works as for B nn . 
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I 1 



10: 
11: 
12: 
13: 

14: 
15: 
16: 
17: 
18: 
19: 
20: 
21: 
22: 
23: 
24: 
25: 
26: 
27: 
28: 
29: 
30: 
31: 
32: 
33: 
34: 
35: 
36: 



hij, af , i,j = 1,2 and M 
Output: B dd given by vectors ri,r2 € M M 
n = [0 : 7? dd /(A/ - 1) : ilf ] 

r 2 (l) = Rf, 2* =0 
for k = 2 : M 

7 * = 2 , 'i( fc ) - 1 

Solve (82) (83) using Prop.|2] = 

if 7i* < ff?iK?/°i 

Compute x\ using )90t 

Compute iEj using j9U 

T (0) _ r r *i^i 
x 12 — l^lJ^i 

Z = 
repeat 

if 



E 21> T21 



si^S) < s 2 (a;§) 



A : 

else 

A : 

end 

Determine step size t 

•^12 ~ [^12 ^^]-JLi 
/ <- I + 1 

12) ~~ ^( x 12 ) 

T22 = * ( a; 12 ) 
Compute X22 using (84) 
else 

72*2 = 
end 

'f 7ai > 7J2 



until 



< € 



: 721- 

: 7| 2 - 



7 2 
else 

7| 
end 

Compute id* 

Compute 

r 2 (fc) = J R 2 



using 1681 
= log 2 (1 



■7j) 



end 



TABLE IV 
Numerical method to compute B dd 
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Fig. 2. Rate regions for gn = g 22 = 1, 312 = 921 = 2, ki = K2 = 0.3 
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VII. Numerical Illustrations 



Fig. 3. Rate regions for g\\ = 922 = 1, 912 = 921 = 2, k\ 



0.85 



In this section, we illustrate how the channel parameters gij , 
Ki, i,j = 1,2 affect the shape of the rate regions. By choosing 
these parameters in a controlled way, instead of randomly 
drawing channel vectors, we can illustrate interesting relations 
between the parameters and the four rate regions. 

In Figs. [2}f4] we illustrate the scenario where the channel 
gains are symmetric with g\\ = 1722 = 1 and 912 = .921 = 2. 
That is, the crosstalk channels are stronger than the direct 
channels. In Fig. [2] we have m = K2 = 0.3, which 
corresponds to a low spatial correlation of amongst the direct 
and crosstalk channels. We see that, even though the crosstalk 
channel gains are high, lZ nn is almost rectangular and all the 
other regions are contained in lZ nn . In this case, there is no 
need of MUD; it costs too much in terms of useful signal 
power to create extra interference in order to enable MUD. 

We illustrate the other extreme case in Fig. [3] Here we 
have Ki = K2 = 0.85, which corresponds to the case where 
the angle between the direct and crosstalk channel vectors is 
small. We see that all the other regions are contained in lZ dd . 
The combination of strong crosstalk channels and high spatial 
correlation, entails that the cost of boosting interference in 
order enable MUD is very small. 




Fig. 4. Rate regions for gn 
k 2 = 0.3 



922 



1. 912 



921 



2, ki = 0.85, 



In Fig.|U we depict the case of K\ = 0.85 and K2 = 0.3, i.e., 
the channels from TXi have high spatial correlation, whereas 
the spatial correlation of the channels from TX2 is low. In this 



12 



Submitted to IEEE TRANS. SIGNAL PROCESS. 



case we have 1Z = lZ nd . The reason is that RX2 experiences 
high interference and has no problem to decode it. On the 
other hand RXi experiences low interference, so it is better to 
treat it as noise. 

For both Fig. H and Fig. H we see that K dd C K dn [j K nd . 
This is something that we frequently observe when the 
channels are i.i.d. Rayleigh and SNR is around dB. The 
explanation is that for lZ dd both links have to sacrifice part of 
the desired signal power in order to enable MUD. 

VIII. Conclusion 

We proposed an efficient method to compute the Pareto 
boundary of the rate region for the two-user MISO IC with 
MUD-capable receivers. The merit of the proposed method, 
compared to the state-of-the-art, is that it avoids the brute-force 
search over all potentially PO beamforming vector pairs. The 
proposed method is efficient in the sense that its complexity 
is constant with respect to the number of transmit antennas. 
For each point on the Pareto Boundary, we solve the quasi- 
concave optimization either by solving a cubic equation or 
running a scalar line search. Using a non-optimized MATLAB 
implementation, we observed that the numerical method to 
compute B nn is as fast as the closed-form method. The 
reason is that, in each iteration of the numerical method in 
Tab. |U we can initiate the iteration with the solution from 
the previous iteration. Therefore, the gradient search requires 
only a few iterations to converge. Finally, we illustrated with 
a few numerical examples, how the spatial correlation of 
the channels affects the shape of the rate region and which 
decoding strategies to use in order to be Pareto optimal. 

It appears unlikely that there is any structure left in the 
problem that can be exploited in order to further improve the 
computational complexity. The fact that the method is fast 
makes it useful to draw illustrations but more importantly as 
a tool for large scale simulations that involve MISO ICs. Our 
method could be used for computing the Pareto boundary of 
the MISO broadcast channel as well. However for that task 
there will also be a need of performing an extra line search to 
find the optimal power allocation. Perhaps, the methodology 
we brought forward here can be applied to other problems too. 
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